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Интеллектуальная система оптимального 
параметрического синтеза источников 
магнитных полей с ферромагнитными 
элементами конструкции 


В работе предложена интеллектуальная система оптимального параметрического синтеза источников 
магнитных полей, содержащих ферромагнитные элементы конструкции как в предположении о постоянстве 
намагниченности полюсных наконечников, Так и для случая учета зависимости намагниченности от 
напряженности магнитного поля. При проведении синтеза использовался разработанный авторами 
гибридный алгоритм оптимизации роем частиц с эволюционным формированием состава популяции, 
позволяющий находить глобальное решение задачи. 


При проектировании источников магнитных полей (ИМП) в подавляющем 
большинстве случаев решается задача анализа распределения магнитного поля, 
причем заданными являются конструктивные параметры магнитной системы (МС). 
Однако больший интерес представляет решение задачи проектирования как задачи 
синтеза, т.е. обратной задачи поиска оптимального размещения в пространстве 
токозадающих катушек и геометрии МС, обеспечивающих априори заданную кон- 
фигурацию магнитного поля в рабочей зоне. Если МС представляется лишь сово- 
купностью катушек, то задача параметрического синтеза обычно сводится к отыска- 
нию значений их оптимальных размеров и взаимного расположения, а также величин 
протекающих в них токов [1]. Задача синтеза значительно усложняется при наличии 
в конструкции МС ферромагнитных элементов (ФЭ), особенно в случае учета их 
нелинейных магнитных свойств. При этом анализ поля даже для одного варианта 
геометрии МС требует существенных затрат машинного времени. Решение же задачи 
синтеза предполагает многократный анализ различных прототипов конструкции. На дан- 
ный момент универсального программного обеспечения (ПО), позволяющего решать 
рассматриваемую задачу оптимального параметрического синтеза ИМП, содер- 
жащих ФЭ, насколько известно авторам, не существует. 

Целью данной работы является создание интеллектуальной системы (ИС) 
оптимального параметрического синтеза ИМП с ФЭ, обеспечивающей требуемую 
топографию поля в рабочем объеме МС. 

Задача параметрического синтеза ставится и решается как задача оптимизации 
некоторой целевой функции, значение которой позволяет оценить степень соответствия 
полученного распределения поля в объеме априори заданному закону. Разработанная ИС 
предполагает представление полюсных наконечников МС в виде наборов ферро- 
магнитных полюсных элементов. При этом в качестве варьируемых параметров 
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используются геометрические размеры катушек, полюсных наконечников, а также 
габариты самого магнитопровода, которые определяют структуру и геометрию всей 
МС в целом, и, как следствие, профиль полюсных наконечников электромагнита. 

В такой постановке решение задачи синтеза МС предполагает многократное 
решение задачи анализа, а эффективность проведения синтеза зависит как от особен- 
ностей применяемого алгоритма оптимизации, так и от используемых методов расчета 
поля. Классические методы оптимизации, традиционно используемые для этих целей, 
ориентированы на определение локального решения задачи поиска экстремума. 
Существенно больший интерес представляет поиск глобального решения. Поэтому 
имеет смысл применение алгоритмов оптимизации, обладающих данными поисковыми 
свойствами. В качестве такого инструмента использовался разработанный авторами 
гибридный алгоритм оптимизации роем частиц с эволюционным формированием 
состава популяции [2], [3], основанный на парадигме роевого интеллекта с использо- 
ванием эволюционных операторов, заимствованных из генетических алгоритмов. Данный 
бионический мультиагентный алгоритм обладает высокой скоростью сходимости, не 
требует вычисления производных оптимизируемой функции, позволяет легко учесть 
ограничения, накладываемые на пространство поиска, и обеспечивает нахождение гло- 
бального оптимума. 

В разработанной ИС параметрического синтеза входными параметрами являются: 
тип МС, количество токозадающих катушек, число полюсных ФЭ, а также ограничения 
на пределы изменений конструктивных параметров МС, магнитные свойства ферро- 
магнетика. Выходными параметрами служат оптимальные значения геометрических 
размеров МС и найденные значения конструктивных параметров полюсных элементов. 

Расчет магнитного поля в рабочем объеме МС предполагает построение математи- 
ческой модели проектируемой системы. В простейшем случае, когда намагниченность 
полюсных наконечников можно считать постоянной и не зависящей от напряжен- 
ности внешнего магнитного поля, что характерно для магнитотвердых материалов, 
обладающих высоким значением коэрцитивной силы, расчет поля, создаваемого 
электромагнитом, может проводиться в соответствии с формулой 


ы 9% Е 
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где Н(О) - напряженность магнитного поля, создаваемого в точке наблюдения О; 
М№ -— количество ферромагнитных полюсных элементов; с; = М.и; — поверхностная 


плотность магнитных зарядов, размещаемых на 7-м полюсном элементе; М = сопя — 
намагниченность полюсных наконечников, п; — внешняя нормаль к поверхности 


полюсного элемента 5;; 7 РО ^ вектор, проведенный из точки истока Р в точку 


наблюдения О. 

При таких допущениях оптимальный параметрический синтез проводился для 
случаев, когда в рабочем объеме МС требуется создание высокооднородного магнит- 
ного поля, а полюсные наконечники представлялись наборами цилиндрических, кони- 
ческих или призматических полюсных элементов [4-8]. Верификация разработанного ПО 
для случая двух цилиндрических полюсных элементов осуществлялась путем сравнения 
полученных при синтезе оптимальных значений их параметров с известными соот- 
ношениями для шим Роуза (0,02%). В случае синтеза аксиально-симметричной МС 
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при использовании цилиндрических полюсных элементов удалось обеспечить 
однородность поля на уровне 20 — 30 ррт (миллионных долей относительного от- 
клонения). Применение в составе полюсов МС конических полюсных элементов 
позволило повысить однородность поля до значения 10 - 20 ррт. А при синтезе С- 
образного электромагнита с призматическими полюсными элементами прямоугольного 
поперечного сечения однородность магнитного поля составила 50 — 130 ррт. 
Достигнутые характеристики распределения поля в рабочем объеме МС вполне 
удовлетворяют требованиям, предъявляемым ко многим техническим устройствам. 
Однако такая идеализация магнитных свойств является приемлемой в ограничен- 
ном количестве приложений. Поэтому математическая модель должна учитывать 
зависимость намагниченности ферромагнетика от напряженности магнитного поля. 
В случае предположения о линейной зависимости намагниченности ферромагнитного 
материала от напряженности поля, которая имеет место для слабых магнитных полей, 
определение распределения поверхностной плотности магнитных зарядов © на ФЭ МС 
возможно с применением аппарата граничных интегральных уравнений (ГИУ) [9]. 
При этом для аксиально-симметричных МС предполагается, что поверхность, ограни- 
чивающая объем, занимаемый ферромагнетиком, разбивается на кольцевые элементы 
в предположении кусочно-постоянного распределения поверхностной плотности 
магнитных зарядов на них, а анализ топографии магнитного поля МС удобно про- 
водить в цилиндрической системе координат. В этом случае распределение повер- 
хностной плотности магнитных зарядов описывается интегральным уравнением 


Фредгольма П рода 
(0) =^.| К(О, Р)о(Р)@Ь +2^поНо(О), ©) 
1 
ав, = _ р 
где ^= па. К(О,Р) - ядро интегрального уравнения; и, Х — магнит- 
и+ + 


ная проницаемость и магнитная восприимчивость ферромагнитного материала; / — 
кривая, задающая контур аксиально-симметричной МС; Но(О) - напряженность 


внешнего поля, создаваемая токозадающими катушками. После нахождения распре- 
деления поверхностной плотности магнитных зарядов о, поле в рабочем объеме МС 
может быть рассчитано по формуле 


ыы в 1 с(Р)РРО ве 
Нч ей 45р +Но(О). (3) 


Решение интегрального уравнения (2) с учетом введенных допущений сводится к 
решению соответствующей системы линейных алгебраических уравнений (СЛАУ) 
большого порядка, в котором число неизвестных достигает десятков тысяч и для 
решения которой применялся метод СМВЕ$ с проецированием на подпространства 
Крылова. Значительного сокращения времени, затрачиваемого на расчет коэффици- 
ентов СЛАУ, удалось достичь путем учета трансляционной симметрии для кольцевых 
элементов разбиения, что кроме выигрыша во времени счета позволяет также 
экономить используемую при проведении вычислений оперативную память компьютера. 
Разработанное на базе ГИУ и гибридного алгоритма оптимизации ПО успешно ис- 
пользовалось для синтеза аксиально-симметричной МС с цилиндрическими полюс- 
ными элементами, предназначенной для генерации высокооднородного магнитного 
поля [9]. При найденных значениях оптимальных параметров полюсных элементов 
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однородность магнитного поля, которую удалось достичь в рабочем объеме, составляет 
около 80 — 90 ррт. Проверка корректности расчета поля при этом проводилась для 
случая бесконечной пластины, намагничиваемой однородным магнитным полем, для 
которой известна аналитическая зависимость между напряженностью внешнего поля 
и величиной поля внутри нее. Погрешность расчета в этом случае не превышала 0,5%. 

И, наконец, учет нелинейной зависимости между намагниченностью ферромагнит- 
ных элементов конструкции и напряженностью магнитного поля, что характерно для 
большинства технических приложений и, особенно, для случаев высокой напря- 
женности полей, можно выполнить с использованием пространственного интеграль- 
ного уравнения [10] 


: 5 1 ы 1 р 
И О АТр + Но(О), (4) 


решаемого совместно с уравнением М =х(Н)Н ‚, описывающим нелинейную зависи- 


мость между величиной намагниченности и напряженностью магнитного поля. В случае, 
когда МС обладает аксиальной симметрией, объем полюсных наконечников и 
магнитопровода Г удобно разбивать на кольцевые элементы и, в предположении 
постоянства намагниченности в пределах каждого элемента разбиения, решение 
уравнения (4) сводить к решению системы нелинейных уравнений большого порядка. 
Верификация разработанного модуля расчета поля в этом случае проводилась также 
для случая бесконечной пластины, которая моделировалась тонким диском, и для 
шара, находящихся в однородном магнитном поле. Погрешность расчета в первом 
случае составляла 0,2%, а во втором — 1,3%. Созданная ИС применялась для 
проектирования коэрцитиметра с частично замкнутой МС [10-12]. В этих работах 
для описания геометрии МС использовалась технология Оп/О{, для ускорения 
процесса синтеза учитывалась трансляционная симметрия элементов разбиения, а 
решение нелинейной системы уравнений большого порядка проводилось с помощью 
метода Ньютона-СМКЕ$. Разработанный в [10] коэрцитиметр позволяет при найден- 
ных значениях конструктивных параметров генерировать магнитное поле с однород- 
ностью на уровне 0,3%. Приведенные примеры свидетельствуют о высокой эф- 
фективности разработанной ИС. Следует заметить, что созданное ПО позволяет 
получать не только поля с высокой степенью однородности, но и поля, имеющие 
более сложную пространственную конфигурацию, например, эта@Н), = сопэй [12]. 


Выводы 


Разработанный программный комплекс позволяет проводить оптимальный 
синтез источников поля с ферромагнитными элементами для различных случаев 
зависимости намагниченности ферромагнетика от напряженности поля и обеспечить 
требуемую конфигурацию магнитного поля в рабочем объеме проектируемых 
электротехнических устройств. 
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О.М. Якимов, В.Я. Гальченко, Д.Л. Остапущенко 

Тнтелектуальна система оптимального параметричного синтезу джерел магнтних пол!в з феро- 
магн гними елементами конструкцИ 

У робот! розглянуто проектування 1нтелектуально! системи оптимального параметричного синтезу 
джерел магнтних пол1в, що вмипують феромагетн! елементи конструкций, у припущенн! посийнот 
намагн!ченост! полюсних наконечник, а також для випадку урахування залежност! намагн1ченост! 
в1д напруженост! магытного поля. При проведенн! синтезу використовувався розроблений авторами 
г1бридний алгоритм оптимзацй роем часток з еволющйним формуванням составу популяций, що 
дозволяе знаходити глобальне ритення задач1. 


А.М. Уамтоу, И. Уа. НёсйепкКо, О.Г. Озаризсйепко 

ПцеШрепЕ Зу$ет ог Орйта]-Рагатейлс Зупе$15 о? Маопейс Ее 4 Зоигсез УУ/НВ Еегготаспейс 
НМетепе о? Сопзбгисвоп 

№ Фе агафе дезлепте оЁ пиеШеепЕ зузет оЁ орйтпа]-рагатлепис зуп®ез15 оЁ тазпейс Не!4 зоигсез, сощашие 
Гегготаэпейс е]етет, ш аззигтарйоп оРреппапеп тазпейтайоп ап4 Юг сазез ассоиийпе оЁ дереп4депсе Бебмееп 
таопейнтайоп ап4 таопейнс ицепзИу, Вауе Бееп сопз1аегед. Пиппе зупез15 ргосез$ Вубиа а1еог т деуеюреа Бу 
Фе аиог$ \’аз 15е4. К аПозуед +0 Нпа г]оБа!| зоаНоп оР Фе ргоШет. 
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